3.72 \(\int \frac{\sqrt{c+d x^2} \left (e+f x^2\right )^{3/2}}{a+b x^2} \, dx\)

Optimal. Leaf size=403 \[ \frac{c^{3/2} \sqrt{e+f x^2} (b e-a f)^2 \Pi \left (1-\frac{b c}{a d};\tan ^{-1}\left (\frac{\sqrt{d} x}{\sqrt{c}}\right )|1-\frac{c f}{d e}\right )}{a b^2 \sqrt{d} e \sqrt{c+d x^2} \sqrt{\frac{c \left (e+f x^2\right )}{e \left (c+d x^2\right )}}}+\frac{f x \sqrt{c+d x^2} (-3 a d f+b c f+4 b d e)}{3 b^2 d \sqrt{e+f x^2}}+\frac{\sqrt{e} \sqrt{f} \sqrt{c+d x^2} (5 b e-3 a f) F\left (\tan ^{-1}\left (\frac{\sqrt{f} x}{\sqrt{e}}\right )|1-\frac{d e}{c f}\right )}{3 b^2 \sqrt{e+f x^2} \sqrt{\frac{e \left (c+d x^2\right )}{c \left (e+f x^2\right )}}}-\frac{\sqrt{e} \sqrt{f} \sqrt{c+d x^2} (-3 a d f+b c f+4 b d e) E\left (\tan ^{-1}\left (\frac{\sqrt{f} x}{\sqrt{e}}\right )|1-\frac{d e}{c f}\right )}{3 b^2 d \sqrt{e+f x^2} \sqrt{\frac{e \left (c+d x^2\right )}{c \left (e+f x^2\right )}}}+\frac{f x \sqrt{c+d x^2} \sqrt{e+f x^2}}{3 b} \]

[Out]

(f*(4*b*d*e + b*c*f - 3*a*d*f)*x*Sqrt[c + d*x^2])/(3*b^2*d*Sqrt[e + f*x^2]) + (f
*x*Sqrt[c + d*x^2]*Sqrt[e + f*x^2])/(3*b) - (Sqrt[e]*Sqrt[f]*(4*b*d*e + b*c*f -
3*a*d*f)*Sqrt[c + d*x^2]*EllipticE[ArcTan[(Sqrt[f]*x)/Sqrt[e]], 1 - (d*e)/(c*f)]
)/(3*b^2*d*Sqrt[(e*(c + d*x^2))/(c*(e + f*x^2))]*Sqrt[e + f*x^2]) + (Sqrt[e]*Sqr
t[f]*(5*b*e - 3*a*f)*Sqrt[c + d*x^2]*EllipticF[ArcTan[(Sqrt[f]*x)/Sqrt[e]], 1 -
(d*e)/(c*f)])/(3*b^2*Sqrt[(e*(c + d*x^2))/(c*(e + f*x^2))]*Sqrt[e + f*x^2]) + (c
^(3/2)*(b*e - a*f)^2*Sqrt[e + f*x^2]*EllipticPi[1 - (b*c)/(a*d), ArcTan[(Sqrt[d]
*x)/Sqrt[c]], 1 - (c*f)/(d*e)])/(a*b^2*Sqrt[d]*e*Sqrt[c + d*x^2]*Sqrt[(c*(e + f*
x^2))/(e*(c + d*x^2))])

_______________________________________________________________________________________

Rubi [A]  time = 0.944809, antiderivative size = 403, normalized size of antiderivative = 1., number of steps used = 7, number of rules used = 7, integrand size = 32, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.219 \[ \frac{c^{3/2} \sqrt{e+f x^2} (b e-a f)^2 \Pi \left (1-\frac{b c}{a d};\tan ^{-1}\left (\frac{\sqrt{d} x}{\sqrt{c}}\right )|1-\frac{c f}{d e}\right )}{a b^2 \sqrt{d} e \sqrt{c+d x^2} \sqrt{\frac{c \left (e+f x^2\right )}{e \left (c+d x^2\right )}}}+\frac{f x \sqrt{c+d x^2} (-3 a d f+b c f+4 b d e)}{3 b^2 d \sqrt{e+f x^2}}+\frac{\sqrt{e} \sqrt{f} \sqrt{c+d x^2} (5 b e-3 a f) F\left (\tan ^{-1}\left (\frac{\sqrt{f} x}{\sqrt{e}}\right )|1-\frac{d e}{c f}\right )}{3 b^2 \sqrt{e+f x^2} \sqrt{\frac{e \left (c+d x^2\right )}{c \left (e+f x^2\right )}}}-\frac{\sqrt{e} \sqrt{f} \sqrt{c+d x^2} (-3 a d f+b c f+4 b d e) E\left (\tan ^{-1}\left (\frac{\sqrt{f} x}{\sqrt{e}}\right )|1-\frac{d e}{c f}\right )}{3 b^2 d \sqrt{e+f x^2} \sqrt{\frac{e \left (c+d x^2\right )}{c \left (e+f x^2\right )}}}+\frac{f x \sqrt{c+d x^2} \sqrt{e+f x^2}}{3 b} \]

Antiderivative was successfully verified.

[In]  Int[(Sqrt[c + d*x^2]*(e + f*x^2)^(3/2))/(a + b*x^2),x]

[Out]

(f*(4*b*d*e + b*c*f - 3*a*d*f)*x*Sqrt[c + d*x^2])/(3*b^2*d*Sqrt[e + f*x^2]) + (f
*x*Sqrt[c + d*x^2]*Sqrt[e + f*x^2])/(3*b) - (Sqrt[e]*Sqrt[f]*(4*b*d*e + b*c*f -
3*a*d*f)*Sqrt[c + d*x^2]*EllipticE[ArcTan[(Sqrt[f]*x)/Sqrt[e]], 1 - (d*e)/(c*f)]
)/(3*b^2*d*Sqrt[(e*(c + d*x^2))/(c*(e + f*x^2))]*Sqrt[e + f*x^2]) + (Sqrt[e]*Sqr
t[f]*(5*b*e - 3*a*f)*Sqrt[c + d*x^2]*EllipticF[ArcTan[(Sqrt[f]*x)/Sqrt[e]], 1 -
(d*e)/(c*f)])/(3*b^2*Sqrt[(e*(c + d*x^2))/(c*(e + f*x^2))]*Sqrt[e + f*x^2]) + (c
^(3/2)*(b*e - a*f)^2*Sqrt[e + f*x^2]*EllipticPi[1 - (b*c)/(a*d), ArcTan[(Sqrt[d]
*x)/Sqrt[c]], 1 - (c*f)/(d*e)])/(a*b^2*Sqrt[d]*e*Sqrt[c + d*x^2]*Sqrt[(c*(e + f*
x^2))/(e*(c + d*x^2))])

_______________________________________________________________________________________

Rubi in Sympy [A]  time = 131.911, size = 360, normalized size = 0.89 \[ \frac{f x \sqrt{c + d x^{2}} \sqrt{e + f x^{2}}}{3 b} - \frac{c^{\frac{3}{2}} f \sqrt{e + f x^{2}} \left (3 a f - 5 b e\right ) F\left (\operatorname{atan}{\left (\frac{\sqrt{d} x}{\sqrt{c}} \right )}\middle | - \frac{c f}{d e} + 1\right )}{3 b^{2} \sqrt{d} e \sqrt{\frac{c \left (e + f x^{2}\right )}{e \left (c + d x^{2}\right )}} \sqrt{c + d x^{2}}} + \frac{\sqrt{e} \sqrt{f} \sqrt{c + d x^{2}} \left (3 a d f - b c f - 4 b d e\right ) E\left (\operatorname{atan}{\left (\frac{\sqrt{f} x}{\sqrt{e}} \right )}\middle | 1 - \frac{d e}{c f}\right )}{3 b^{2} d \sqrt{\frac{e \left (c + d x^{2}\right )}{c \left (e + f x^{2}\right )}} \sqrt{e + f x^{2}}} - \frac{f x \sqrt{c + d x^{2}} \left (3 a d f - b c f - 4 b d e\right )}{3 b^{2} d \sqrt{e + f x^{2}}} + \frac{c^{\frac{3}{2}} \sqrt{e + f x^{2}} \left (a f - b e\right )^{2} \Pi \left (1 - \frac{b c}{a d}; \operatorname{atan}{\left (\frac{\sqrt{d} x}{\sqrt{c}} \right )}\middle | - \frac{c f}{d e} + 1\right )}{a b^{2} \sqrt{d} e \sqrt{\frac{c \left (e + f x^{2}\right )}{e \left (c + d x^{2}\right )}} \sqrt{c + d x^{2}}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  rubi_integrate((d*x**2+c)**(1/2)*(f*x**2+e)**(3/2)/(b*x**2+a),x)

[Out]

f*x*sqrt(c + d*x**2)*sqrt(e + f*x**2)/(3*b) - c**(3/2)*f*sqrt(e + f*x**2)*(3*a*f
 - 5*b*e)*elliptic_f(atan(sqrt(d)*x/sqrt(c)), -c*f/(d*e) + 1)/(3*b**2*sqrt(d)*e*
sqrt(c*(e + f*x**2)/(e*(c + d*x**2)))*sqrt(c + d*x**2)) + sqrt(e)*sqrt(f)*sqrt(c
 + d*x**2)*(3*a*d*f - b*c*f - 4*b*d*e)*elliptic_e(atan(sqrt(f)*x/sqrt(e)), 1 - d
*e/(c*f))/(3*b**2*d*sqrt(e*(c + d*x**2)/(c*(e + f*x**2)))*sqrt(e + f*x**2)) - f*
x*sqrt(c + d*x**2)*(3*a*d*f - b*c*f - 4*b*d*e)/(3*b**2*d*sqrt(e + f*x**2)) + c**
(3/2)*sqrt(e + f*x**2)*(a*f - b*e)**2*elliptic_pi(1 - b*c/(a*d), atan(sqrt(d)*x/
sqrt(c)), -c*f/(d*e) + 1)/(a*b**2*sqrt(d)*e*sqrt(c*(e + f*x**2)/(e*(c + d*x**2))
)*sqrt(c + d*x**2))

_______________________________________________________________________________________

Mathematica [C]  time = 2.92614, size = 739, normalized size = 1.83 \[ \frac{3 i a^3 d f^2 \sqrt{\frac{d x^2}{c}+1} \sqrt{\frac{f x^2}{e}+1} \Pi \left (\frac{b c}{a d};i \sinh ^{-1}\left (\sqrt{\frac{d}{c}} x\right )|\frac{c f}{d e}\right )-i a \sqrt{\frac{d x^2}{c}+1} \sqrt{\frac{f x^2}{e}+1} \left (3 a^2 d f^2-3 a b f (c f+d e)+b^2 e (4 c f-d e)\right ) F\left (i \sinh ^{-1}\left (\sqrt{\frac{d}{c}} x\right )|\frac{c f}{d e}\right )-3 i a^2 b c f^2 \sqrt{\frac{d x^2}{c}+1} \sqrt{\frac{f x^2}{e}+1} \Pi \left (\frac{b c}{a d};i \sinh ^{-1}\left (\sqrt{\frac{d}{c}} x\right )|\frac{c f}{d e}\right )-6 i a^2 b d e f \sqrt{\frac{d x^2}{c}+1} \sqrt{\frac{f x^2}{e}+1} \Pi \left (\frac{b c}{a d};i \sinh ^{-1}\left (\sqrt{\frac{d}{c}} x\right )|\frac{c f}{d e}\right )-3 i b^3 c e^2 \sqrt{\frac{d x^2}{c}+1} \sqrt{\frac{f x^2}{e}+1} \Pi \left (\frac{b c}{a d};i \sinh ^{-1}\left (\sqrt{\frac{d}{c}} x\right )|\frac{c f}{d e}\right )+3 i a b^2 d e^2 \sqrt{\frac{d x^2}{c}+1} \sqrt{\frac{f x^2}{e}+1} \Pi \left (\frac{b c}{a d};i \sinh ^{-1}\left (\sqrt{\frac{d}{c}} x\right )|\frac{c f}{d e}\right )+a b^2 d e f x^3 \sqrt{\frac{d}{c}}+6 i a b^2 c e f \sqrt{\frac{d x^2}{c}+1} \sqrt{\frac{f x^2}{e}+1} \Pi \left (\frac{b c}{a d};i \sinh ^{-1}\left (\sqrt{\frac{d}{c}} x\right )|\frac{c f}{d e}\right )+a b^2 c e f x \sqrt{\frac{d}{c}}+a b^2 d f^2 x^5 \sqrt{\frac{d}{c}}+a b^2 c f^2 x^3 \sqrt{\frac{d}{c}}-i a b e \sqrt{\frac{d x^2}{c}+1} \sqrt{\frac{f x^2}{e}+1} (-3 a d f+b c f+4 b d e) E\left (i \sinh ^{-1}\left (\sqrt{\frac{d}{c}} x\right )|\frac{c f}{d e}\right )}{3 a b^3 \sqrt{\frac{d}{c}} \sqrt{c+d x^2} \sqrt{e+f x^2}} \]

Antiderivative was successfully verified.

[In]  Integrate[(Sqrt[c + d*x^2]*(e + f*x^2)^(3/2))/(a + b*x^2),x]

[Out]

(a*b^2*c*Sqrt[d/c]*e*f*x + a*b^2*d*Sqrt[d/c]*e*f*x^3 + a*b^2*c*Sqrt[d/c]*f^2*x^3
 + a*b^2*d*Sqrt[d/c]*f^2*x^5 - I*a*b*e*(4*b*d*e + b*c*f - 3*a*d*f)*Sqrt[1 + (d*x
^2)/c]*Sqrt[1 + (f*x^2)/e]*EllipticE[I*ArcSinh[Sqrt[d/c]*x], (c*f)/(d*e)] - I*a*
(3*a^2*d*f^2 - 3*a*b*f*(d*e + c*f) + b^2*e*(-(d*e) + 4*c*f))*Sqrt[1 + (d*x^2)/c]
*Sqrt[1 + (f*x^2)/e]*EllipticF[I*ArcSinh[Sqrt[d/c]*x], (c*f)/(d*e)] - (3*I)*b^3*
c*e^2*Sqrt[1 + (d*x^2)/c]*Sqrt[1 + (f*x^2)/e]*EllipticPi[(b*c)/(a*d), I*ArcSinh[
Sqrt[d/c]*x], (c*f)/(d*e)] + (3*I)*a*b^2*d*e^2*Sqrt[1 + (d*x^2)/c]*Sqrt[1 + (f*x
^2)/e]*EllipticPi[(b*c)/(a*d), I*ArcSinh[Sqrt[d/c]*x], (c*f)/(d*e)] + (6*I)*a*b^
2*c*e*f*Sqrt[1 + (d*x^2)/c]*Sqrt[1 + (f*x^2)/e]*EllipticPi[(b*c)/(a*d), I*ArcSin
h[Sqrt[d/c]*x], (c*f)/(d*e)] - (6*I)*a^2*b*d*e*f*Sqrt[1 + (d*x^2)/c]*Sqrt[1 + (f
*x^2)/e]*EllipticPi[(b*c)/(a*d), I*ArcSinh[Sqrt[d/c]*x], (c*f)/(d*e)] - (3*I)*a^
2*b*c*f^2*Sqrt[1 + (d*x^2)/c]*Sqrt[1 + (f*x^2)/e]*EllipticPi[(b*c)/(a*d), I*ArcS
inh[Sqrt[d/c]*x], (c*f)/(d*e)] + (3*I)*a^3*d*f^2*Sqrt[1 + (d*x^2)/c]*Sqrt[1 + (f
*x^2)/e]*EllipticPi[(b*c)/(a*d), I*ArcSinh[Sqrt[d/c]*x], (c*f)/(d*e)])/(3*a*b^3*
Sqrt[d/c]*Sqrt[c + d*x^2]*Sqrt[e + f*x^2])

_______________________________________________________________________________________

Maple [B]  time = 0.027, size = 1028, normalized size = 2.6 \[ \text{result too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  int((d*x^2+c)^(1/2)*(f*x^2+e)^(3/2)/(b*x^2+a),x)

[Out]

1/3*(d*x^2+c)^(1/2)*(f*x^2+e)^(1/2)*((-d/c)^(1/2)*x^5*a*b^2*d*f^2+(-d/c)^(1/2)*x
^3*a*b^2*c*f^2+(-d/c)^(1/2)*x^3*a*b^2*d*e*f+3*((d*x^2+c)/c)^(1/2)*((f*x^2+e)/e)^
(1/2)*EllipticF(x*(-d/c)^(1/2),(c*f/d/e)^(1/2))*a^3*d*f^2-3*((d*x^2+c)/c)^(1/2)*
((f*x^2+e)/e)^(1/2)*EllipticF(x*(-d/c)^(1/2),(c*f/d/e)^(1/2))*a^2*b*c*f^2-3*((d*
x^2+c)/c)^(1/2)*((f*x^2+e)/e)^(1/2)*EllipticF(x*(-d/c)^(1/2),(c*f/d/e)^(1/2))*a^
2*b*d*e*f+4*((d*x^2+c)/c)^(1/2)*((f*x^2+e)/e)^(1/2)*EllipticF(x*(-d/c)^(1/2),(c*
f/d/e)^(1/2))*a*b^2*c*e*f-((d*x^2+c)/c)^(1/2)*((f*x^2+e)/e)^(1/2)*EllipticF(x*(-
d/c)^(1/2),(c*f/d/e)^(1/2))*a*b^2*d*e^2-3*((d*x^2+c)/c)^(1/2)*((f*x^2+e)/e)^(1/2
)*EllipticE(x*(-d/c)^(1/2),(c*f/d/e)^(1/2))*a^2*b*d*e*f+((d*x^2+c)/c)^(1/2)*((f*
x^2+e)/e)^(1/2)*EllipticE(x*(-d/c)^(1/2),(c*f/d/e)^(1/2))*a*b^2*c*e*f+4*((d*x^2+
c)/c)^(1/2)*((f*x^2+e)/e)^(1/2)*EllipticE(x*(-d/c)^(1/2),(c*f/d/e)^(1/2))*a*b^2*
d*e^2-3*((d*x^2+c)/c)^(1/2)*((f*x^2+e)/e)^(1/2)*EllipticPi(x*(-d/c)^(1/2),b*c/a/
d,(-f/e)^(1/2)/(-d/c)^(1/2))*a^3*d*f^2+3*((d*x^2+c)/c)^(1/2)*((f*x^2+e)/e)^(1/2)
*EllipticPi(x*(-d/c)^(1/2),b*c/a/d,(-f/e)^(1/2)/(-d/c)^(1/2))*a^2*b*c*f^2+6*((d*
x^2+c)/c)^(1/2)*((f*x^2+e)/e)^(1/2)*EllipticPi(x*(-d/c)^(1/2),b*c/a/d,(-f/e)^(1/
2)/(-d/c)^(1/2))*a^2*b*d*e*f-6*((d*x^2+c)/c)^(1/2)*((f*x^2+e)/e)^(1/2)*EllipticP
i(x*(-d/c)^(1/2),b*c/a/d,(-f/e)^(1/2)/(-d/c)^(1/2))*a*b^2*c*e*f-3*((d*x^2+c)/c)^
(1/2)*((f*x^2+e)/e)^(1/2)*EllipticPi(x*(-d/c)^(1/2),b*c/a/d,(-f/e)^(1/2)/(-d/c)^
(1/2))*a*b^2*d*e^2+3*((d*x^2+c)/c)^(1/2)*((f*x^2+e)/e)^(1/2)*EllipticPi(x*(-d/c)
^(1/2),b*c/a/d,(-f/e)^(1/2)/(-d/c)^(1/2))*b^3*c*e^2+(-d/c)^(1/2)*x*a*b^2*c*e*f)/
(d*f*x^4+c*f*x^2+d*e*x^2+c*e)/b^3/(-d/c)^(1/2)/a

_______________________________________________________________________________________

Maxima [F]  time = 0., size = 0, normalized size = 0. \[ \int \frac{\sqrt{d x^{2} + c}{\left (f x^{2} + e\right )}^{\frac{3}{2}}}{b x^{2} + a}\,{d x} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  integrate(sqrt(d*x^2 + c)*(f*x^2 + e)^(3/2)/(b*x^2 + a),x, algorithm="maxima")

[Out]

integrate(sqrt(d*x^2 + c)*(f*x^2 + e)^(3/2)/(b*x^2 + a), x)

_______________________________________________________________________________________

Fricas [F(-1)]  time = 0., size = 0, normalized size = 0. \[ \text{Timed out} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  integrate(sqrt(d*x^2 + c)*(f*x^2 + e)^(3/2)/(b*x^2 + a),x, algorithm="fricas")

[Out]

Timed out

_______________________________________________________________________________________

Sympy [F]  time = 0., size = 0, normalized size = 0. \[ \int \frac{\sqrt{c + d x^{2}} \left (e + f x^{2}\right )^{\frac{3}{2}}}{a + b x^{2}}\, dx \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  integrate((d*x**2+c)**(1/2)*(f*x**2+e)**(3/2)/(b*x**2+a),x)

[Out]

Integral(sqrt(c + d*x**2)*(e + f*x**2)**(3/2)/(a + b*x**2), x)

_______________________________________________________________________________________

GIAC/XCAS [F]  time = 0., size = 0, normalized size = 0. \[ \int \frac{\sqrt{d x^{2} + c}{\left (f x^{2} + e\right )}^{\frac{3}{2}}}{b x^{2} + a}\,{d x} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  integrate(sqrt(d*x^2 + c)*(f*x^2 + e)^(3/2)/(b*x^2 + a),x, algorithm="giac")

[Out]

integrate(sqrt(d*x^2 + c)*(f*x^2 + e)^(3/2)/(b*x^2 + a), x)